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Abstract— We  develop  a  general  theorem  for  the  force  on  singularities  in  an  elastic  two-dimensional 
medium  in  terms  of  the  residues  of  the  stress  at  the  singularities.  This  theorem  is  then  applied  to  a  sharp 
crack  screened  by  a  cloud  of  dislocations.  A  total  fracture  criterion  can  then  be  derived  in  principle  by 
specifying  the  local  cleavage  condition  at  the  crack  and  the  lattice  resistance  of  the  dislocations.  The 
COD  of  the  crack  is  shown  to  be  given  by  the  total  screening  Burgers  vector  of  the  dislocation  cloud 
and  the  wake  of  a  moving  crack  is  discussed  in  terms  of  the  resistance  to  moving  the  screening  cloud. 
Finally,  limitations  of  the  model  are  discussed  in  terms  of  the  geometrical  effects  on  the  cleavage  plane 
caused  by  blunting  the  crack. 

Khnmt — Nous  diveloppons  un  theorems  general  concemant  la  force  sur  des  singularity  dans  un 
milieu  ilastique  bidimensionnel  4  partir  des  residus  de  la  contrainte  sur  ces  singularity.  On  applique  ce 
theoreme  4  une  fissure  etroite  6crantee  par  un  nuage  de  dislocations.  En  principe,  on  peut  obtenir 
un  entire  de  rupture  totale  en  pricisant  la  condition  locale  de  clivage  sur  la  fissure  et  la  resistance 
du  resea u  sur  les  dislocations.  On  montre  que  le  nuage  de  la  fissure  est  defini  par  le  vecteur  de 
Burgers  d'ecrantage  total  de  ce  nuage  de  dislocations  et  Ton  discute  la  mise  en  mouvement  d'une  fissure 
en  fonction  de  la  resistance  que  Ton  doit  surmonter  pour  deplacer  cette  atmosphere.  Enfin.  on  discutc 
les  limitations  du  modele  en  termes  des  edicts  giometriques  sur  le  plan  de  clivage.  provoques  par  une 
fissure  qui  s’emousse. 

Zusammenfanuag  Wir  formutieren  ein  allgcmeines  Theorem  fiir  die  Kraft,  die  auf  Singularitaten  in 
einem  elastischen  zweidimensionalen  Medium  wirkt,  auf  der  Basis  der  Residucn  der  Spannung  an  den 
Singularitaten.  Dieses  Theorem  wird  auf  einen  scharfen  RiB,  der  von  einer  Versetzungswolke 
abgeschirmt  ist.  angewendet.  Ein  Bruchkriterium  kann  im  Prinzip  abgelcitet  werden,  wenn  die  lokalen 
Spaltungsbedingungen  am  RiB  und  der  Gleitwiderstand  fiir  die  Versetzungen  angegeben  wird.  Die 
kritische  RiBdffnungs-verschiebung  its  bestimmt  durch  den  Gesamtburgersvektor  der 
Versetzungswolke.  Der  Verlauf  eines  sich  bewegenden  Risses  wird  anhand  des  Widerstandes  diskutiert, 
den  die  Abschirmwolke  bei  ihrer  Bewegung  uberwinden  muB.  Zuletzt  werden  die  Grenzcn  des  Modelles 
an  den  geometrischen  Effekten  auf  der  Spaltungsebene.  die  durch  Abstumpfung  des  Risses  entstehen. 
aufgezeigt. 


1.  INTRODUCTION 

In  this  paper  we  shall  develop  a  general  theorem  for 
the  force  on  singularities  in  an  elastic  two-dimen¬ 
sional  medium  in  terms  of  the  residues  of  the  stress  at 
the  singularity  (in  the  sense  of  complex  function 
theory),  and  apply  this  theorem  to  a  sharp  crack 
shielded  by  a  distribution  of  dislocations  in  anti-plane 
strain  (mode  III  fracture).  From  this  application,  a 
dislocation  basis  for  fracture  mechanics  can  be  devel¬ 
oped  which  permits  more  satisfactory  approaches  to 
fracture  criteria  than  is  possible  through  continuum 
fracture  mechanics  alone. 

The  theory  is  based  on  the  properties  of  cracks  and 
dislocations  described  as  singularities  in  an  elastic 
field,  and  thus  addresses  only  the  behavior  of  sharp 
cracks  with  their  associated  shielding  clouds  of  dislo¬ 
cations.  In  the  paper  we  shall  nevertheless  discuss 
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how  a  significant  crack  opening  displacement  (COD) 
can  be  inferred  for  such  a  crack  and  how  sharp  cracks 
governed  by  cleavage  processes  at  their  tips  are 
related  to  the  fully  ductile  variety. 

Our  motivation  for  developing  such  a  theory  based 
on  discrete  elastic  singularities  is  that  in  continuum 
theory  there  is  no  satisfactory  way  to  develop  a  frac¬ 
ture  criterion  based  on  the  microstructurc  of  the 
material.  Further,  in  the  ./-integral  theory,  the 
medium  is  required  to  possess  an  energy  density  func¬ 
tion  which  is  single  valued  at  all  points.  This  require¬ 
ment  makes  it  impossible  to  consider  a  moving  crack 
rigorously,  because  as  the  crack  moves,  those  portions 
of  the  medium  sufficiently  close  to  the  crack  are  first 
loaded  up  as  the  crack  approaches,  undergo  plastic 
deformation  and  finally  unload  as  the  crack  passes. 
Since  by  definition,  plastic  deformation  is  highly  hys- 
teretic.  such  a  loading  cycle  will  violate  the  assump¬ 
tion  of  single  valuedness,  and  make  the  /-integral  in- 
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applicable.  This  limitation  is  well  known  in  the  frac¬ 
ture  mechanics  community  [1]  and  its  implications 
have  been  intensively  studied,  because  the  concept  of 
specifying  the  critical  onset  of  fracture  in  terms  of  Jlc 
requires  taking  account  of  the  effects  of  unloading  [2]. 

In  a  second  example,  if  the  maximum  yield  strength 
at  the  tip  of  the  crack  is  finite,  as  it  must  be  in  any 
real  material.  Rice  [1]  has  shown  that  the  force  on  the 
bare  crack  is  effectively  zero,  because  all  the  work 
done  during  crack  growth  is  entirely  absorbed  by  the 
plastic  deformation,  and  that  a  fracture  criterion  can¬ 
not  be  constructed  by  analogy  to  the  Griffith  energy 
balance  concepts  of  linear  elastic  brittle  fracture.  This 
is  an  intuitively  unfortunate  result,  because  then  the 
crack  opening  processes  cannot  be  separated  physi¬ 
cally  or  mathematically  from  the  plastic  work  pro¬ 
cesses  in  the  ‘process  zone’. 

On  the  other  hand,  if  one  could  build  a  theory 
which  described  the  deformation  aspects  of  fracture  in 
terms  of  the  behavior  of  the  dislocations  which  in  fact 
constitute  the  plastic  zone,  one  might  hope  to  circum¬ 
vent  these  difficulties  which  are  inherent  in  the  con¬ 
tinuum  approximation.  Further,  because  a  fracture 
mechanics  based  on  dislocation  concepts  would  be 
more  fundamentally  grounded,  one  might  hope  to 
build  a  framework  within  which  to  address  more 
satisfactorily  such  questions  as  the  materials  par¬ 
ameters  controlling  the  fracture  criterion,  and  the 
time  and  temperature  properties  of  fracture,  at  least 
qualitatively.  We  are  guided  by  the  hope  that  by 
focussing  on  the  dislocation  aspects  of  fracture,  one 
may  be  able  to  develop  an  alternative  picture  which 
can  usefully  complement  the  usual  continuum  ap¬ 
proximation. 

The  beginnings  of  a  fracture  mechanics  based  on 
dislocation  theory  have  already  been  made.  Rice  [3] 
has  derived  an  expression  for  the  force  on  a  dislo¬ 
cation  near  a  crack,  and  for  the  change  of  the  stress 
intensity  factor  of  the  crack  induced  by  the  dislo¬ 
cation.  Thomson  [4],  Weertman  [5]  and  Thomson 
and  Fuller  [6]  have  introduced  the  idea  of  a  dislo¬ 
cation  shielded  crack,  with  an  elastic  enclave  near  the 
tip  and  Hart  [7]  has  developed  a  continuum  dislo¬ 
cation  screening  theory.  Majumdar  and  Burns  [8] 
have  derived  stress  distributions  and  discussed  inter¬ 
action  forces  for  a  variety  of  crack-dislocation  con¬ 
figurations.  Finally.  Chang  and  Ohr[9]  have  dis¬ 
cussed  the  dislocation  free  zone  surrounding  a  crack 
tip. 

2.  THE  FORCE  ON  AN  ELASTIC 
SINGULARITY 

In  this  section  we  show  how  the  expression  for  the 
force  on  an  elastic  singularity  in  terms  of  a  contour 
integral  of  Eshelby's  energy-momentum  tensor  assumes 
a  particularly  simple  form  in  the  complex-variable 
notation  possible  for  two-dimensional  linear-clastic 
fields.  This  result  is  given  here  for  the  case  of 
general  anisotropy,  and  can  be  used  to  develop  simple 
expressions  for  energy  release  rates  associated  with 


motion  of  defects  of  all  kinds.  Budiansky  and 
Rice  [10]  have  derived  corresponding  expressions  for 
the  forces  on  simple  cracks  for  the  isotropic  plane  and 
anti-plane  strain  cases,  but  they  have  not  put  their 
results  into  the  form  of  the  residues  at  the  singulari¬ 
ties  as  we  shall  do,  nor  have  they  included  dislo¬ 
cations  in  their  treatment.  In  subsequent  sections  of 
this  paper  we  shall  use  the  extremely  simple  form  for 
the  anti-plane  strain  case. 

For  a  field  of  displacements  u,  with  associated 
stress  Ojj,  dependent  only  on  x,  and  x2,  Eshelby  [II] 
has  derived  the  expression 

9i  =  <j>  Pt/ijds  (1) 

for  the  force  on  a  singularity  or  singularities  con¬ 
tained  within  the  contour  S  in  the  x,x2  plane,  where 
ny  is  the  outward  unit  normal  to  S,  and  the  energy- 
momentum  tensor  Ptj  is  given,  for  linear  elasticity,  by 

P ,j  =  i cruUt.Ay  —  OkjUk.i-  (2) 

The  usual  suffix  notations  for  summations  and  deriva¬ 
tives  are  assumed.  These  forces  may  be  caused  by 
stresses  exerted  on  the  external  surfaces  of  the 
medium,  or  alternatively  by  other  singularities  such  as 
dislocations,  cracks,  etc.,  which  happen  to  lie  outside 
S.  When  applied  to  a  crack,  g  is  the  same  as  Irwin's 
crack  extension  force.  StJ  is  the  stress,  and  u,  the  dis¬ 
placement.  In  the  usual  way,  this  theorem  is  applied 
to  elastic  cracks  by  noting  that  P,/nJ  +  «}*)  =  0, 
where  n'j  and  n'/  refer  to  normal  vectors  on  opposite 
sides  of  the  mathematical  slit  representing  the  cleav¬ 
age  surface.  Thus  jP^fiy  ds  =  0.  when  integrated  over 
the  two  cleavage  surfaces,  provided  the  fracture  sur¬ 
faces  are  free  of  applied  forces.  Thus,  a  contour  for  a 
slit  crack  which  ends  at  opposite  points  of  the  cleav¬ 
age  surface  is  equivalent  to  one  which  follows  the 
cleavage  surface  to  the  crack  tip  and  back  (see  Fig.  1). 

Eshelby  et  at.  [12]  first  gave  general  complex-vari¬ 
able  expressions  for  displacement  and  stress  fields  in 
equilibrium  under  anisotropic  elasticity. 

In  the  variant  notation  of  Stroh[13]  these  ex¬ 
pressions  take  the  form 


U,  =  2Rej  i  At,f,(i,) J 

(3a) 

(  3  1 

®i,  =  2Rej-  I  P.W,'<z,)j 

(3b) 

f  3  1 

=  2Re  )j 

(3c) 

where  /,(:>) are  analytic  functions  of  z,  =  x,  +  P,x2. 
/,(z»)  s  d/,(z,)/dz,  and  the  complex  constants  P„ 
A„  and  Li,  are  functions  of  the  elastic  constants. 
After  manipulations  involving  relations  between  A  „ 
and  Li,  not  given  here,  the  required  components  of 
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Fig.  1.  A  contour  around  a  slit  crack.  Since  the  integral  for 
g  over  the  cleavage  surface  (ADC)  contained  within  S  is 
zero  when  no  forces  are  exerted  on  that  portion  of  the 
cleavage  surface,  S  may  be  considered  to  consist  only  of  the 
portion,  ABC.  S  may  enclose  other  singularities  beside  the 
crack,  and  still  others  may  lie  external  to  S. 


done  during  crack  growth  is  entirely  absorbed  by  the 
Pu  are  found  to  be 


Pu  =  -Pu  -  2 Re  j  i  P,Lta/l*.t/^,)32J 

(4a) 

P12  =  2Rej-  £  LtaA*[/;<z,)]2  J 

(4b) 

P2,  =  2Re{j£  P2Lb/tb[/;(z,)]2J. 

(4c) 

Substitution  into  (1)  yields 

g,  =  2Re|  X  F,.Lte/4„ 

x  (  ±  j  U'.MV  <»*) 

}  ,5' 

where 


P„  ^1  Pu  m  P, 

the  ±  sign  is  the  sign  of  Im(P,),  and  S,  are  the 
contours  in  the  Z,  planes  corresponding  to  S. 
Provided  S  lies  entirely  in  elastic  material  in 
equilibrium  the  (/i)2  are  analytic  along  S.  and  (5) 
can  be  expressed 

ft  -  4*lm  J  £  ±  PhL+A* 

*  IRes[/,;<2.)]2J  (6) 

where  IRes  denotes  the  sum  of  the  residues  of  the 
squared  stress  function  at  all  its  poles  inside  S,. 

(n  the  limit  of  isotropy,  the  plane  and  anti-plane 
strain  problems  decouple:  anti-plane  solutions  are 
given  in  terms  of  just  one  stress  function,  /3.  We  have 

P3— »i;  Lt3— *  — *  -Sk)/2fi 

where  p  is  the  shear  modulus,  equations  (3b)  and  (3c) 
become 

(7) 


it,  where  z  =  x  +  iy,  and  equations  (S)  and  (6)  become 

0  =  0i  -  »0j  =  2^:  (8) 

=  ^£Res(<72).  (9) 

As  an  elementary  example,  for  a  semi-infinite  slit 
crack  on  x2  =  0,  x,  <  0,  loaded  in  mode  III, 

e (z)  =  K/yj2itz  (10) 

where  K  is  the  stress  intensity  factor.  In  this  case,  n1 
is  analytic  even  on  the  negative  real  axis  where  a  is 
discontinuous,  and  (9)  gives  for  the  crack  extension 
force 

0,  =K2/ '2n.  (11) 

Similarly,  for  an  isolated  screw  dislocation  under  a 
remote  uniform  stress, 

a  =  a,  +  iA>/2n:  (12) 

and  the  force  is  given  by 

g  =  <T|b.  (13) 

Although  these  simple  results  are  familiar,  the  real 

advantage  of  the  present  formulation  becomes  more 
apparent  when  collections  of  defects,  such  as  cracks 
with  accompanying  dislocations,  are  considered.  The 
additional  terms  in  the  complex  stress  functions 
necessary  to  satisfy  the  combined  boundary  con¬ 
ditions  can  be  readily  identified  with  interaction  and 
image  forces.  Strain  energy  release  rates  for  cooper¬ 
ative  motion  of  defects  can  be  readily  evaluated  by 
summing  residues  at  the  associated  singularities.  The 
celebrated  path  independence  of  the  ./-integral  (  =g, ) 
is  clearly  identified  with  the  path  independence  of 
contour  integrals  of  meromorphic  functions  provided 
the  poles  enclosed  remain  constant  (Cauchy  residue 
theorem). 


3.  STRESS  AND  FORCE  CALCULATIONS 
FOR  SCREENED  CRACKS 

We  now  wish  to  apply  the  above  theorem  to  the 
case  of  mode  III  fracture  with  contained  yield,  treat¬ 
ing  the  screening  dislocations  as  discrete  entities. 

We  postulate  a  crack  which  in  the  complex  plane 
x  +  iy  occupies  the  negative  real  axis,  and  at  first 
consider  just  one  dislocation,  at  the  point  f,  Fig.  2. 
With  full  anisotropy.  Sinclair  and  Hirth  [14]  have 

j.t 


Fig.  2.  Schematic  of  coordinate  system  in  r  -  x  +  iy 
plane.  A  dislocation  at  C  exerts  stresses  on  the  cleavage 
surface  which  must  be  cancelled  by  a  distribution  of 
sources  at  on  the  cleavage  surface. 


ffu  +  iffj,  -/i(z)  a  tr(z) 
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crack  and  its  screen  by  considering  contours,  S,  which 
enclose  these  and  only  these  singularities.  First, 
however,  we  write  the  familiar  expression  for  a  simple 
bare  crack  without  any  dislocation  screen  present. 
That  is,  the  total  stress  is  given  by  equation  (10).  Then 
for  any  contour  surrounding  the  origin 

K> 

0  bare  cr.ch  *  -  (17) 


Fig.  3.  The  crack  plus  dislocation  screen.  S,  encloses  only 
the  crack  tip,  and  leads  to  g, ;  S2  surrounds  a  single  dislo¬ 
cation  in  thus  screen,  and  leads  to  gd\  and  S3  surrounds  the 

entire  configuration,  and  leads  to  g,„Al. 

solved  the  corresponding  problem  for  a  finite  length 
crack,  and  this  can  be  taken  to  the  limit  of  infinite 
crack  length  after  a  suitable  coordinate  shift  to  center 
one  tip  at  the  origin.  Alternatively,  Majumdar  and 
Burns  [8]  have  solved  the  present  case  directly  in 
anti-plane  strain.  By  either  route,  we  obtain  for  the 
complex  stress  function. 

4 njz  ly'r  -  v/f  J z  +  vf J 


With  the  screen  present.  Fig  3,  and  evaluating  equa¬ 
tion  (9)  with  (IS)  for  the  residue  at  z  =  0  correspond¬ 
ing  to  a  contour  surrounding  the  origin,  but  enclosing 
no  dislocations,  we  obtain  the  force  on  the  core  crack 
by  itself  in  the  presence  of  the  screen 


This  result  states  that  the  crack  extension  force  is  the 
same  as  the  Irwin-Orowan  result,  except  that  the  core 
crack  is  shielded  from  the  external  stress  field  by  its 
dislocation  screen.  It  is  suggestive  to  write  this  as 

9,  =  ~  *2-.  r  =  K-K„ 


=  /*bf_J _ L_ 

4rr  (z  -  (  z  -  ( 

In  the  expanded  second  form,  the  first  term  in  the 
braces  is  that  for  a  dislocation  without  crack:  the 
second  is  an  image  term,  the  third  represents  the  fact 
that  any  external  stress  applied  to  the  crack  generates 
a  crack-like  singularity  at  the  tip  and  the  fourth  is  the 
image  term  for  the  third. 

We  can  now  obtain  the  stress  at  any  point  z  due  to 
the  combination  of  (1)  a  slit  core  crack  at  the  origia 
and  (2)  screening  dislocations  distributed  at  discrete 
points,  Zj,  by  adding  the  stresses  of  the  individual  con¬ 
tributions 


We  emphasize  that  gc  is  the  local  force  to  extend  the 
core  crack,  independent  of  the  motion  of  the  dislo¬ 
cation  screen.  The  shielding  effect  predicted  by  equa¬ 
tion  (19)  has  been  discussed  in  various  forms  by  all 
the  previous  authors  listed  in  Refs  [3-8]. 

The  force  on  one  of  the  dislocations  of  the  screen  is 
obtained  from  equations  (15)  and  (9)  by  identifying 
the  residue  at  the  position  of  one  of  the  dislocations, 
say  that  at  C-  Then  for  a  contour  surrounding  C  and 
no  other  singular  points,  we  have 
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(20) 

Equation  (15)  reduces  to  a  simpler  form  if  the  distri¬ 
bution  is  symmetric  about  the  x-axis.  i.e.  if  for  every 
dislocation  at  £/  there  is  a  matching  one  at  £y  with  the 
same  b. 


g  ^  +  y  fa  _ 

V'2trr  )  2*  \J  :  z  - 


(16) 


This  is  the  result  given  in  integral  form  by  Hart  [7], 
Having  a  rigorous  expression  for  the  stress  at  any 
point  we  can  now  apply  (15)  to  obtain  the  force  on 
any  element  or  collection  of  elements  of  the  core 


Fig.  4.  Schematic  stress-strain  law.  o  vs  «.  with  r„  the 
critical  yield  stress. 


l 
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The  summation  sign  is  primed  to  remind  us  that 
the  sum  is  not  taken  over  the  dislocation  at  (■  This 
expression  can  also  be  written  in  terms  of  * 


(Mi) 
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(21) 


The  second  form  is  obviously  an  expansion  about  the 
near  region  of  the  crack  tip.  where  t  is  dominant, 
while  the  first.  (20),  is  an  expansion  about  the  outer 
edge  of  the  screen,  where  K  is  dominant. 

By  a  direct  summation  over  the  residues  at  the 
crack  and  the  dislocation  positions,  using  a  contour 
surrounding  the  entire  configuration,  we  can  also 
write  an  expression  for  the  total  force  on  the  entire 
crack-dislocation  configuration.  This  we  will  term  the 
‘dressed  crack’ 


g,M.i  =  *  I  Res(tr2(o)  +  <r2(y) 
^  J 

=  KI/2ft 


(22) 


where  g,„„  is  the  force  to  move  the  entire  configur¬ 
ation  rigidly  through  the  medium.  Note  it  is  the  same 
as  one  would  calculate  if  the  dressed  crack  configur¬ 
ation  was  an  effective  bare  crack  with  a  simple 
K-singularity  at  the  origin  as  in  equation  (17). 

It  is  instructive  to  compare  these  results  with  those 
of  a  continuum  deformation  theory.  In  continuum 
theory,  the  stress-strain  law  is  a  nonlinear  function  at 
stresses  above  r0  (Fig.  4),  but  the  law  is  single  valued 
without  hysteresis,  as  explained  in  the  introduction. 
Under  these  circumstances,  the  ./-integral  is  applic¬ 
able,  and  is  contour  independent.  The  ./-integral  is 
defined  by 

J  s  g.  (23) 

where  gx  is  defined  in  equation  (8).  In  the  elastic 
region,  the  equation 

a  =  Kv  2n:  (24) 

holds  and  J  can  be  calculated  on  a  contour  where  (24) 
is  valid.  Thus 

J  =  =  K2  2  g.  (25) 

We  note  that  even  if  the  stress-strain  law  is  perfectly 
plastic  for  2«  >  t0  t‘.  the  7-integral  is  finite  for  an 
infinitesimal  contour  around  the  crack  tip.  because 
the  strain  is  singular  at  the  tip. 

The  connection  between  these  continuum  results 
and  ours  is  that,  unlike  J,  our  g  is  not  invariant  as  the 
contour  is  shrunk  through  the  deformation  region 
containing  dislocations  (see  Fig.  5).  The  reason  is  that 
as  the  contour  is  shrunk,  small  circles  are  left  behind 


(•) 


(b) 


PLASTIC  ZONE 


Fig.  $.  In  (a),  the  deformation  is  represented  by  actual  dis¬ 
locations.  Thus  as  S,  is  shrunk  to  S2  it  leaves  behind  a 
wake  of  enclosed  singularities  if  g  is  to  remain  invariant. 
Sj,  since  it  encloses  a  different  set  of  poles  than  S,  and  S2, 
does  not  have  the  same  value.  In  (b).  the  deformation  in  the 
shaded  region  is  modeled  as  a  nonlinear  stress-strain  re¬ 
lation.  and  g  is  an  invariant:  giS,)  =  gfS2). 


surrounding  each  of  the  dislocations  outside  the  main 
contour.  When  we  draw  a  contour  leaving  out  all  or 
some  portion  of  these  singularities,  the  g  thus  calcu¬ 
lated  will  be  different  from  that  using  a  contour  in  the 
elastic  region  surrounding  the  entire  dressed  crack. 
Contour  integrations  thus  remain  invariant  when  the 
contours  are  shrunk  (or  deformed  in  some  other  way) 
only  in  so  far  as  they  do  not  pass  through  one  of  the 
singularities  of  the  configuration.  This  breakdown  of 
the  path  independence  theorem  for  the  screened 
cracks  will  be  seen  to  have  important  consequences. 


4.  IMPLICATIONS  FOR  FRACTURE 
MECHANICS 

The  mathematical  theorems  in  the  previous  two 
sections  will  now  be  applied  to  some  of  the  central 
problems  of  fracture  mechanics.  From  our  standpoint, 
the  ‘solution’  of  the  fracture  problem  is  obtained 
when  one  is  able  to  find  a  self-consistent  distribution 
of  dislocations  (perhaps  time  dependent)  such  that 
there  is  zero  net  force  (including  other  material 
induced  forces  specified  below)  on  all  elements  of  the 
total  configuration.  This  is  an  enormous  task,  of 
course,  and  is  not  to  be  attempted  directly  except  in 
some  very  special  cases.  Otherwise,  it  will  be  necess¬ 
ary  to  develop  suitable  approximate  approaches. 
However,  there  are  some  important  rigorous  results 
and  insights  which  can  be  obtained  without  carrying 
through  such  a  program  in  detail. 
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4.1.  Governing  equations  for  crack  equilibrium 

As  explained  above,  a  static  configuration  is 
achieved  only  when  each  singularity  is  in  local  static 
force  equilibrium.  The  forces  derived  in  [18-21]  rep¬ 
resent  crack  extension  forces  from  external  stresses, 
and  dislocation-dislocation  and  dislocation-crack  in¬ 
teraction  forces.  However,  the  lattice  itself  exerts  ad¬ 
ditional  forces  on  the  crack  and  dislocations  due  to 
surface  tension,  interactions  with  chemical  inhomogen¬ 
eities,  and  lattice  resistance  from  lattice  trapping  and 
Peierls  forces.  When  these  additional  terms  are  in¬ 
cluded,  we  can  write  the  basic  governing  equations  for 
static  equilibrium  of  the  configuration  as 

*3/2g  =  2y  +  t 

0a(O  =  rt.  (26) 

The  first  equation  represents  local  equilibrium  for  the 
crack  in  terms  of  the  surface  tension  or  energy,  y  and 
the  lattice  trapping  barriers,  simply  written  as  t.  In 
essence,  it  is  simply  an  expression  of  the  Griffith  re¬ 
lation  for  the  core  crack,  because  in  true  static  equi¬ 
librium  t  =  0  and  the  t  terms  are  only  important  for 
slow  crack  growth.  Although  it  has  not  been  ob¬ 
served,  except  in  some  special  cases,  it  should  be  kept 
in  mind  that  materials  and  circumstances  exist  for 
which  trapping  may  be  sufficiently  large  that  t  must 
be  included  explicitly  in  the  Griffith  equilibrium  con¬ 
dition.  For  the  dislocations  we  have  simply  written 
the  resistance  terms  as  r4  is  intended  to  include 
the  complex  set  of  material  phenomena  such  as  im¬ 
purity  pinning,  grain  boundary  effects  and  Peierls 
barriers  which  must  be  juggled  by  the  theorist  in 
modeling  all  dislocation  mobility  problems.  The  first 
equation  in  (26)  is  of  particular  interest,  because  it 
expresses  a  necessary  fracture  criterion  for  sharp 
cracks  which  is  beyond  anything  derivable  from  con¬ 
tinuum  plasticity  theories.  In  continuum  plastic 
theories,  the  asymptotic  stress  dependence  at  the 
crack  tip  is  a  power  law  different  from  - 1/2  and  it 
has  not  been  possible  to  link  these  theories  up  to  the 
cohesive  forces  at  the  crack  tip  which  hold  the  atoms 
of  the  crack  tip  together  [I],  On  the  other  hand,  (26) 
follows  in  our  theory  from  the  fact  that  the  dislo¬ 
cations  are  discrete,  and  the  core  crack  has  a  K  field 
at  its  tip.  This  result  can  only  be  self-consistent  if  the 
density  of  dislocations  near  the  tip  is  low  enough  that 
dislocation  cores  and  crack  cores  do  not  overlap.  It 
must  also  be  true  that  the  crack  cleavage  surface  is 
not  blunted  by  the  action  of  the  dislocations,  but  we 
will  discuss  that  separately  in  a  later  section. 

Equation  (26)  also  predicts  that  an  overalt  fracture 
criterion  will  feature  the  intrinsic  surface  energy,  y.  in 
an  important  way.  It  is  often  said  that  y  is  such  a 
small  portion  of  the  'effective  surface  energy'  for  a 
moving  crack  that  the  overall  criterion  will  be  insensi¬ 
tive  to  it.  However,  y  is  one  of  the  important  force 
balancing  parameters  driving  the  whole  fracture  pro¬ 
cess.  because  it  specifies  the  essential  core  crack  equi¬ 
librium  condition  and  crack  tip  stress  intensity  factor. 


.nr  to  which  the  rest  of  the  dislocation  configuration 
must  conform.  Thus,  a  change  in  y  will  generate  first 
order  changes  in  the  overall  equilibrium  criterion. 
Attempts  to  calculate  the  critical  K,c  for  crack 
growth,  by  approximate  means  confirm  this  predic¬ 
tion  [4-6], 

Equation  (26)  refers  to  a  static  crack.  However, 
brittle  cracks  in  some  materials  are  known  to  be  able 
to  sustain  stable  crack  growth  over  a  range  of  K 
values.  In  such  a  case,  specifying  the  value  of  K  speci¬ 
fies  a  velocity  for  the  core  crack.  Time  dependence 
and  rate  effects  must  then  be  introduced  into  the  re¬ 
sponse  of  the  dislocation  screen.  Hart  [7]  has  cons- 
structed  such  a  quasi-static  theory  on  the  basis  of  a 
continuum  constitutive  creep  law  for  the  dislocation 
screen.  His  work  is  limited  to  dislocation  density  dis¬ 
tributions  which  lead  to  asymptotic  1  /Jz  stresses  at 
the  crack  tip.  Our  results,  however,  show  that  his 
approach  has  a  larger  generality  provided  the  discrete 
character  of  the  dislocation  distribution  near  the 
crack  tip  and  the  core  crack  K  field  are  treated  con¬ 
sistently. 


4.2.  Energy  balance  criteria  and  moving  cracks 

The  equilibrium  conditions  (26)  are  part  of  a  pre¬ 
scription  for  finding  the  detailed  distributions  of  the 
dislocations  in  a  complete  solution  of  the  fracture 
problem.  However,  there  are  also  integral  relations 
which  are  derivable  from  total  energy  conservation  of 
a  system  in  which  an  external  stress  drives  a  crack 
through  a  material  in  steady  state.  Without  asking 
how  the  state  is  attained,  we  assume  the  crack  moves 
quasi-statically  through  the  medium  carrying  its 
screen  with  it.  but  also  leaving  a  wake  of  plastically 
worked  material  behind. 

To  describe  this  situation  adequately  requires  us  to 
distinguish  between  the  geometric  and  nongeometric 
dislocation  density.  In  any  real  material,  during  plas¬ 
tic  deformation,  dislocation  generation  is  a  very 
chaotic  process  in  which  dislocations  of  all  allowed 
Burgers  vectors  are  formed.  Thus  there  is  a  large 
degree  of  cancellation  of  long-range  stress.  To  be 
quantitative,  we  suppose  in  mode  III  that  we  define 
the  local  geometric  Burgers  vector  by  summing  the 
Burgers  vectors  of  the  actual  dislocations  over  an  area 
which  contains  a  large  number  of  dislocations,  and 
divide  by  the  area  of  the  element.  Since  all  by  are  in 
the  X)  direction,  we  have 

I*; 

B  -  =  B(r).  (28) 

Over  cells  of  the  medium  containing  many  dislo¬ 
cations,  B(r)  will  be  a  reasonably  well  defined  and 
continuous  function.  B(r)  is  the  geometric  density,  and 
gives  rise  to  the  long-range  stress  which  we  actually 
wish  to  use  in  all  sums  like  equation  (IS)  to  calculate 
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the  shielding  of  the  crack.  On  the  other  hand 


(29) 


where  N  is  the  actual  number  of  dislocations  thread¬ 
ing  the  area  in  question.  In  general  the  actual  work  of 
fracture  can  be  expected  to  be  from  one  to  several 
orders  of  magnitude  larger  than  the  work  to  form  the 
geometric  density. 

In  these  terms,  as  a  crack  progresses  through  a 
material,  in  the  plastic  zone  dislocations  are  generated 
with  the  necessary  local  geometric  density,  B(z),  to 
shield  the  crack.  In  addition,  in  order  for  the  crack 
and  its  geometric  screening  cloud  to  move  uniformly, 
generation,  motion  and  a  large  degree  of  annihilation 
of  the  real  dislocations  will  be  required.  In  order  to 
simulate  this  complex  system,  we  assume  that  the 
screening  geometrical  dislocation  density,  B(z),  is 
known  and  is  stationary  in  a  reference  system  rigidly 
fixed  to  the  uniformly  moving  crack,  Fig.  6.  In  ad¬ 
dition,  because  of  the  irreversibility  of  the  deforma¬ 
tion,  a  wake  is  generated.  This  wake  is  composed  of 
tangled  dislocations  whose  long-range  stress  is  zero, 
many  of  whom  ultimately  annihilate,  and  which  we 
will  represent  as  a  distribution  of  dislocation  dipoles 
in  the  crystal.  The  wake  contains  a  large  amount  of 
stored  and  thermal  energy,  however,  and  since  it  is 
produced  by  the  deformation  held,  it  exerts  an  effec¬ 
tive  drag  force  on  the  geometric  screening  dis¬ 
locations  against  their  motion  through  the  crystal. 
The  drag  force  on  the  bare  crack  is  y  and  the  drag 
force  on  the  geometric  dislocations  is  given  by  the 
rate  of  formation  of  stored  dislocation  and  thermal 
energy  in  the  wake.  If  the  wake  contains  only  dipoles, 
then  it  exerts  no  elastic  force  on  the  crack  and.  of 
course,  to  a  first  approximation,  the  dislocation 
dipoles  are  themselves  force  free  in  the  (slowly  vary¬ 
ing)  stress  field  of  the  wake. 

Thus  the  total  driving  force  on  crack  plus  screen  is 
K1/ 2  p  and  the  force  balance  becomes 

K1 

—  =  2y  +  R  +  2  T  (30) 

2g 

2 rSx  is  the  stored  dislocation  dipole  energy  in  a 
section  Sx  wide  perpendicular  to  the  cleavage  plane, 
RSx  is  the  energy  dissipated  into  heat  and  y  is  the 
intrinsic  surface  energy.  This  equation  is,  of  course, 
simply  an  expression  of  energy  conservation,  but  it 
cannot  be  obtained  by  straightforward  application  of 
the  ./-integral.  We  noted  earlier  that  >'.  even  though 
small,  is  an  important  term  in  an  e^  ation  like  (30). 
which  means  that  f  and  R  depend  on  y.  because  the 
size  of  the  plastic  zone  depends  upon  y. 

We  note  here  that  there  is  no  reason  to  expect  B(z) 
to  be  strictly  zero  in  the  wake,  but  this  point  is 
explored  in  the  next  section. 


4.3  Burgers  vector  conservation  and  crack  blunting 

Dislocation  sources  in  two  dimensions  create  dislo¬ 
cations  in  pairs  of  dislocations  with  opposite  sign. 
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Fig.  6.  A  uniformly  growing  crack  generates  a  screening 
distribution  in  the  active  deformation  region  surrounding 
the  core  crack  tip  and  leaves  behind  a  debris  of  plastically 
worked  material  in  which  the  long-range  stress  is  zero.  In 
(b)  the  geometrical  dislocations  are  assumed  to  move 
rigidly  with  the  crack,  and  leave  behind  a  wake  of  dislo¬ 
cation  dipoles  which  become  the  plastically  deformed 
wake.  If  the  wake  is  composed  purely  of  dipoles,  the  crack 
faces  in  the  wake  region  are  parallel  and  the  separation  is 
the  COD. 


Hence,  the  total  dislocation  distribution  will  possess 
zero  net  Burgers  vector  and  a  local  excess  of  one  sign 
must  be  balanced  elsewhere  by  an  excess  of  opposite 
sign,  unless  one  class  of  dislocations  is  absorbed  sys¬ 
tematically  by  the  cleavage  surfaces  of  the  crack. 
From  equation  (18),  the  dislocations  which  act  to 
shield  the  crack  from  the  external  stress  are  repelled 
from  the  crack  by  the  fC -field  of  the  crack  at  distances 
beyond  the  influence  of  the  short  range  image  terms. 
Their  partners,  by  equation  (18).  must  be  anti-shield¬ 
ing  dislocations  which  are  attracted  by  the  crack 
K  -field.  Thus  for  the  dislocation  screen  to  have  a  net 
overall  shielding  effect  on  the  crack,  the  anti-shielding 
dislocations  must  be  absorbed  by  the  crack  surface 
(see  Fig.  7). 

Since  anti-shielding  dislocations  which  are  ab¬ 
sorbed  by  the  cleavage  surface  create  slip  steps  of  one 
predominant  sign  on  that  surface,  they  contribute  the 
plastic  component  of  the  crack  opening  displacement. 
By  conservation  of  Burgers  vector,  then,  the  total 
Burgers  vector  in  the  screen  is  equal  to  the  crack 
opening  displacement 


COD  =  £<b,),  (31) 

l 

(this  result  is  limited  to  two  dimensions,  because  in 
three  dimensions,  shielding  can  occur  without  COD 
contributions  when  the  slip  planes  are  skew  to  the 
crack  line.  Thus  (31)  is  only  an  approximate  law  in 
real  materials.) 
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Fig.  7.  Conservation  of  Burgers  vector  at  a  crack.  Initially, 
the  total  number  of  shielding  and  anti-shielding  dislo¬ 
cations  is  equal.  The  anti-shielding  dislocations  are 
attracted  to  the  crack,  and  tend  to  be  annihilated  in  the 
cleavage  surface.  The  COD  is  equal  to  the  total  number  of 
anti-shielding  dislocations  annihilated. 

Although  (31)  has  been  derived  assuming  that  the 
dislocations  are  created  in  pairs  at  internal  sources, 
the  same  result  follows  if  the  dislocations  are  created 
at  the  crack  tip.  because  when  a  dislocation  is  created, 
the  crack  is  blunted  by  one  Burgers  vector. 

Equation  (31)  provides  the  fundamental  connection 
between  the  crack  opening  displacement  and  the 
overall  toughness,  because  a  high  screening  'charge'  of 
total  Burgers  vector  must  be  roughly  associated  with 
higher  toughness.  The  precise  relation,  however, 
depends  upon  the  details  of  the  distribution  of  the 
screen,  i.e.  far  away  dislocations  do  not  screen  as  well 
as  close  in  ones,  etc. 

We  may  make  use  of  the  continuum  solution  for 
the  mode  III  fracture  to  calculate  the  predicted  COD- 
toughness  relation  in  that  case.  From  Ref.  [4],  we 
take  the  approximate  relation 


vector  Zbj,  is  then  given  by 

4>«/2  j»2Xrcoi0 

COD  =  2  dfl  b,(r,  0)r  dr 

Jo  Jo 

=  SR^O  = -g/io-  (34) 

n 

In  deriving  (34)  we  have  used  the  result  K2  =  2nxlRp 
where  Rp  is  the  plastic  zone  radius.  Because  of  the 
approximation  of  placing  the  crack  tip  at  the  edge  of 
the  plastic  circle,  (34)  is  only  valid  for  the  limit  where 
nsO.  This  equation,  however,  is  in  good  agreement 
with  the  BCS  result  as  quoted  in  equation  (218)  of 
Ref.  [16],  and  with  the  various  dimensional  argu¬ 
ments  connecting  J  with  COD  [16].  In  spite  of  the 
relation  (34),  however,  it  is  clear  that  the  rough  equiv¬ 
alence  between  COD  and  g/x0  is  model  dependent 
and  for  extreme  dislocation  distributions  (34)  may  be 
poorly  obeyed. 

We  mentioned  in  the  previous  section  that  the 
wake  may  contain  a  non-zero  Bfz).  One  could  expect 
this  to  be  the  case  often,  because  as  the  crack  passes, 
some  portion  of  the  geometric  component  will  surely 
be  trapped  by  the  wake  dipole  component.  These  dis¬ 
locations  will,  of  course,  be  attracted  to  the  open 
cleavage  surface,  but  only  those  very  near  this  surface 
will  actually  be  able  to  push  their  way  through  the 
work-hardened  region  of  the  wake  to  annihilate  on 
the  open  surface.  The  moving  crack  in  this  case  will 
generate  net  Burgers  vector  on  a  continuous  basis. 
See  Fig.  9.  Thus,  as  the  crack  moves  a  distance  6x,  it 
leaves  behind  a  net  Burgers  vector  given  by 

B«5x  =  f  '  [B(z)]y  dy  <5x  (35) 


<3« 


(32) 


€j„  is  the  elastic  portion  of  the  strain.  The  coordinate 
system  is  depicted  in  Fig.  (8).  The  tip  of  the  crack  is 
placed  at  the  edge  of  the  plastic  zore.  as  in  perfect 
plasticity,  whereas  in  a  work  hardening  material  with 
work  hardening  exponent,  n,  the  tip  of  the  crack 
extends  a  small  distance  into  the  plastic  zone.  See 
Ref.  [4],  e  j,  is  the  elastic  component  of  the  strain,  r0 
and  c0  are  the  critical  shear  stress  and  strain,  respect¬ 
ively.  The  density  of  geometric  Burgers  vector  is  given 
by 


b ,(r,0)  =  2(V  x  <), 

2c  o  /K2  cos  fl  Y’' 141,1 
(I  +  nK  l  ittor  / 

2c  o  /  RpY"'1  *■" 

(l+nVV*/ 


(33) 


(This  expression  for  density  of  Burgers  vector  corrects 
equations  (8)  and  (9)  of  Ref.  [4]).  The  total  Burgers 


where  Rp  is  the  width  of  the  wake.  B  is  the  Burgers 
vector  contained  in  the  wake  per  unit  length  of  wake 
n’ong  one  cleavage  surface.  By  the  arguments  of  this 
section,  the  net  Burgers  vector  is  reflected  in  ad¬ 
ditional  crack  opening,  so  that  if  B(z)  in  the  wake 
does  not  depend  on  x  in  the  wake,  the  crack  shape 
becomes  a  simple  wedge,  with  opening  angle  given  by 
x, 

tan  a  =  B.  (36) 

If  the  wake  retains  some  geometric  component,  it  will 
also  exert  elastic  forces  on  the  crack  and  on  the 
mobile  screening  dislocations  whose  distribution 


Fig.  8.  In  the  continuum  solution,  the  strain  is  a  function 
of  r.  If  work-hardening  is  present,  the  circle  overlaps  the 
crack  tip  slightly  for  small  n.  but  we  approximate  the  result 
by  putting  the  crack  tip  on  the  radius  of  the  circle. 
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Fig.  9.  When  the  wake  retains  some  net  Burgers  vector,  the 
crack  develops  a  wedge  shape  with  an  angle,  a.  related  to 
the  net  geometric  dislocation  density. 


moves  with  the  crack.  If  we  assume  as  in  (35)  that  Rp 
is  the  radius  of  the  'mobile'  screen  and  likewise  the 
width  of  the  wake,  then  the  geometric  dislocations 
Bdx  in  the  wake  at  distance  x  from  the  crack  tip 
where  x  »  Rp  will  actually  appear  to  the  crack  as  a 
dipole,  because  of  the  image  term  of  B  in  the  cleavage 
surface.  Thus,  except  for  the  near  region  of  the  wake 
where  .x  ==  Rr  the  elastic  interaction  of  the  wake  with 
both  the  screening  charge  and  the  crack  is  negligible. 

4.4  Mixed  continuum-discrete  approaches 

As  noted  earlier  the  complete  self-consistent  sol¬ 
ution  of  the  fracture  problem  is  not  feasible  in  prac¬ 
tice,  and  recourse  must  be  taken  to  approximate 
methods.  Clearly,  deep  in  the  plastic  zone  where  a 
continuous  dislocation  density  function  can  be 
defined,  the  continuum  plastic  approach  is  appro¬ 
priate.  The  usual  constitutive  law  between  stress  and 
plastic  strain  takes  the  place  of  the  governing  equilib¬ 
rium  equations  for  the  dislocations.  The  difficulty 
comes  in  the  matter  of  how  to  cut  off  the  continuum 
solution  near  the  crack  tip,  and  what  kind  of  bound¬ 
ary  condition  to  set  on  the  inner  continuum  bound¬ 
ary.  Two  papers  have  dealt  with  this  question  in  dif¬ 
ferent  ways  [4. 5].  Both  have  postulated  an  inner  elas¬ 
tic  enclave  which  models  the  fact  that  the  crack  at  its 
tip  sees  a  A -field,  and  have  used  the  Griffith  relation 
at  the  tip.  (26).  to  obtain  the  fracture  criterion. 

Two  difficulties  have  emerged  from  this  work.  The 
first  relates  to  the  appropriate  boundary  conditions  to 
specify  on  the  boundary  of  the  elastic  enclave.  This 
problem  is  not  straightforward,  because  the  enclave 
boundary  is  the  locus  of  the  cut-off  for  the  plastic 
solution  (the  dislocation  density  function  loses  its 
meaning,  because  it  degenerates  into  sparsely  spaced 
(5-functions),  but  the  elastic  portion  of  the  solution 
remains  valid  down  to  atomic  dimensions.  Since  elas¬ 
tic  and  plastic  displacements  are  intrinsically  tied  to 
one  another  in  a  continuum  plasticity  theory,  a 
proper  procedure  for  separating  the  one  from  the 
other  must  be  fashioned  on  physical  grounds.  So 
far  [6].  a  fully  satisfactory  method  for  doing  this  has 
not  been  fashioned. 


In  these  theories,  however,  the  size  of  the  elastic 
enclave  is  an  important  parameter,  and  this  is  the 
second  difficulty  referred  to  above.  Little  guidance  has 
been  given  for  making  estimates  of  the  enclave  size.  In 
the  terms  of  this  paper,  the  radius  of  the  elastic 
enclave  is  equal  to  the  scale  of  dislocation  hetero¬ 
geneity  at  the  crack  tip;  i.e.  given  by  the  limiting  den¬ 
sity  of  dislocations,  limiting  spacing  of  glide  planes, 
etc.  It  is  precisely  on  this  point  where  the  discreteness 
of  the  dislocation  description  is  crucial,  but  also 
where  our  knowledge  is  most  inadequate,  because  so 
little  is  known  about  the  characteristics  of  very  dense 
dislocation  distributions  in  regions  of  combined  high 
stress  and  high  stress  gradient. 

In  addition  to  the  expected  role  to  be  played  by  the 
dislocation-dislocation  interactions  and  source  den¬ 
sity  and  placement  in  determining  the  enclave  radius, 
it  is  possible  that  the  image  vs  K -field  competition  in 
equation  (21)  plays  a  crucial  role.  Assume  that  a  dislo¬ 
cation  is  very  near  the  crack,  and  interacts  with  it 
more  strongly  than  with  all  the  other  dislocations. 
Then  the  first  term  and  the  first  brace  in  (21)  are  the 
dominant  terms.  We  note  that  the  two  terms  in  ques¬ 
tion  are  of  opposite  sign,  and  that  one  dominates  at 
close  distances  (the  image  term  in  the  brace)  and  the 
other  at  greater  distance  (the  K  term).  Thus  the  force 
on  a  shielding  dislocation  (repulsive  at  large  distance) 
shows  a  maximum  at  a  few  atomic  distances  from  the 
tip  which  turns  out  be  be  quite  broad.  That  is,  the 
repulsive  term  is  significantly  softened  out  to  the 
order  of  several  tens  of  atomic  distances.  This  is  far 
enough  to  be  a  serious  contender  in  determining  the 
radius  of  the  elastic  enclave.  Physically,  what  happens 
is  that  at  short  distances,  the  crack  tip  loses  its  hard 
core  of  repulsion  against  the  ‘pile  up'  of  shielding 
dislocations  and  this  sets  a  limit  to  the  maximum 
stress  in  the  ‘pile  up'  at  the  crack  tip.  Theoretical 
estimates  have  been  made  of  the  extent  of  the  maxi¬ 
mum  in  Ref.  [3]  and  Majumdar  and  Bums  [15] 
believe  it  has  been  observed  in  experiments  in  both 
ionic  crystals  and  metals.  In  our  opinion,  further 
detailed  studies  of  an  experimental  nature  are  badly 
needed  on  this  point. 
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